epi draft 



Excitations of Ellipsoid Packings near Jamming 



Z. Zeravcic^'^ (a)^ Xu2'3, A. J. Liu^, S. R. Nagel^ and W. van Saarloos^ 

^ Instituut - Lorentz, Universiteit Leiden, Posthus 9506, 2300 RA Leiden, The Netherlands 

James Franck Institute and Department of Physics, University of Chicago, Chicago IL 60637 USA 

(N3 

Department of Physics and Astronomy, University of Pennsylvania, Philadelphia PA 19104 USA 

^ ■ 

PACS 63 . 50 . -X - Vibrational states in disordered systems 
PACS 45 . 70 . -n - Granular systems 
PACS 63 . 50 . Lm - Glasses and amorphous solids 
PACS 61 . 43 . - j - Disordered solids 



O 



o 
u 

(N 
> 
00 

in 
o 

On 
O 



Abstract. - We study the vibrational modes of three-dimensional jammed packings of soft ellip- 
soids of revolution as a function of particle aspect ratio e and packing fraction. At the jamming 
transition for ellipsoids, as distinct from the idealized case using spheres where £ = 1, there are 
many unconstrained and non-trivial rotational degrees of freedom. These constitute a set of zero- 
frequency modes that are gradually mobilized into a new rotational band as |e — 1| increases. Quite 
surprisingly, as this new band is separated from zero frequency by a gap, and lies below the onset 
frequency for translational vibrations, u* , the presence of these new degrees of freedom leaves 
unaltered the basic scenario that the translational spectrum is determined only by the average 
contact number. Indeed, uj* depends solely on coordination as it does for compressed packings of 
spheres. We also discuss the regime of large je — 1|, where the two bands merge. 
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Introduction. — There is no doubt that increasing 
ihc pressure on a floppy assembly of particles can create 
a rigid material. Less obvious is the assertion that such 
matter caused by jamming is generically different from 
ordinary solids. At zero temperature, frietionless, ideal 
spheres jam when the number of inter-partiele contacts 
is precisely that necessary to constrain all motion; this 
leads to an excess of decidedly unusual low-frequency vi- 
brational modes. However, it has been argued that frie- 
tionless spheres, with no significant rotational degrees of 
freedom, represent a singular situation that is too ideal- 
ized to reflect the true nature of the marginally-jammed 
state. In particular, jammed ellipsoids can have dramat- 
ically fewer contacts than needed to inhibit all motion so 
that there will be many zero- frequency vibrations. Thus 
the distribution of collective excitations could be funda- 
mentally altered by particle shape. Contrary to this ex- 
pectation, we find the spherical jamming transition, like 
many singular points, controls a broader class of behavior 
but in an unusual non-trivial way. The new degrees of free- 
dom introduced by ellipsoids create a rotational band but 
leave unaltered the basic scenario that the translational 
spectrum is determined solely by the average number of 
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contacts per particle. 

We study these issues in the context of the spectral 
density of vibrational excitations, D{w). In a three- 
dimensional solid, the low-frequency spectrum should fol- 
low the Debye law Di^uj) ~ dictated by the elastic 
modes. This is one of the most robust generic behaviors 
in all of materials science. However, this law breaks down 
in a spectacular fashion for the case of a rigid solid formed 
from the jamming of spheres interacting via finite-ranged 
repulsions [1~7]. The onset of jamming in such systems has 
features of a first-order transition, with a discontinuity in 
the number of interacting neighbors per particle [1,8], as 
well as features of a second-order transition, with power- 
law scaling and diverging length scales [1-14]. Just above 
the zero-temperature transition, D{uj) is approximately 
constant down to zero frequency [1,2] implying the ex- 
istence of a new class of low-frequency excitations that 
arise because the solid is on the verge of instability [3-5]. 
The Maxwell criterion for rigidity [15] proposes that the 
average number of interacting neighbors per particle, Z , 
should be high enough to constrain all relevant degrees of 
freedom in the sample: Z > Zc- For frietionless spheres, 
the critical coordination number, Zc = 6, coincides [1] 
with the value found at the jamming threshold packing 
fraction, (pc- At packing fractions (j) > (f>c, Z exceeds Zc 
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and consequently the plateau in the density of states per- 
sists only down to a frequency lo* that depends solely on 
SZ = (Z — Zc) [2-5]. The apparent gap emerging in the 
spectrum between lo = and lo = to* contains ordinary 
elastic plane-waves described by Debye theory. Is the new 
physics of the excess modes robust for jamming transitions 
generally [16] or is it applicable only to this idealized sit- 
uation of spheres? 

It was succinctly demonstrated [17] that in one sense 
spheres represent a singular situation and therefore may 
be a poor starting point for describing the generic prop- 
erties of jammed solids. The introduction of even a small 
distortion to a sphere introduces many new degrees of 
freedom that need to be constrained for complete sta- 
bility. While a sphere has only three relevant (transla- 
tional) degrees of freedom, a spheroid (an ellipsoid of rev- 
olution with one symmetry axis) requires two additional 
coordinates (two Euler angles) to specify its orientation. 
Maxwell's counting argument for the rigidity of spheroid 
packings would necessitate an average coordination num- 
ber Zc ~ 10. Clearly a discontinuous increase in den- 
sity would be needed if the introduction of an arbitrarily 
small distortion of the sphere required the average num- 
ber of contacts per particle to jump discontinuously from 
6 to 10. The rapid increase of the coordination number 
Z with distortion and, in particular, of the packing frac- 
tion [17-22] has garnered much attention ("M&M's pack 
more efficiently than spheres" [17,18,23]). Nevertheless, 
at the jamming threshold Z increases smoothly — not 
discontinuously — from Z = 6 as spheres deform into el- 
lipsoids, so that for small distortion Z is below Zc = 10 in 
apparent violation of the Maxwell criterion. While there 
are exactly the minimum number of contacts needed for 
mechanical stability at the jamming transition of spheres, 
there are fewer than the minimum number needed for el- 
lipsoids. There must therefore be unconstrained degrees 
of freedom [20] so that the solid is not stable (to quadratic 
order) to some excitations. Here we investigate how rota- 
tions introduce new non-zero-frequency excitations. Re- 
markably we find that these modes do not destroy the 
picture developed for spheres but instead can be naturally 
incorporated into this scenario. 

The Maxwell stability argument and the occur- 
rence of zero modes. — It is clearly somewhat odd to 
think of Zc as jumping discontinuously as soon as there is a 
minute distortion of the particles from spherical symmetry, 
so it is useful to briefly reconsider the Maxwell's argument 
for the stability/instability threshold [15]. The Maxwell 
criterion is based on global counting arguments for the 
minimum number of contacts needed to maintain force 
balance on all particles. Since contacts are shared by two 
particles, the minimum number of contacts necessary to 
clamp all particles which experience forces, Zc, is accord- 
ing to this argument twice the number of degrees of free- 
dom of a particle. As the apparent discontinuity in Zc sim- 
ply arises from our decision whether or not to include the 



rotational degrees of freedom in the counting, it is more 
intuitive to restore continuity by thinking of each sphere 
as having 6 degrees of freedom: three for translations and 
three for arbitrary rotations so that Zc = 12. The rota- 
tions of individual (frictionless) spheres do not contribute 
in any way to the stability of the packing and are thus 
simply the (Zc — Z)/2 = 3 zero-frequency modes per par- 
ticle that are trivially localized on each particle. A natu- 
ral scenario is that these innocuous zero-frequency modes 
progressively become mobilized into finite-frequency exci- 
tations with increasing distortion of the spheres into el- 
lipsoids. There are clearly two important values of the 
coordination number, Z = 6 and Zc = 12. The important 
issue taken up here is the question: which of these con- 
trols the spectrum of excitations for the generalized case 
of ellipsoids? 

From the above point of view, just above the jamming 
threshold, the fact that Z is below Zc should manifest itself 
in the presence of (Zc ~ Z)/2 normal modes of vibration 
per particle with zero frequency — see Fig. la. We study 
here the nature of these modes and the question of how 
they become mobilized into finite-frequency excitations so 
as to find out whether this process changes the jamming 
scenario of frictionless spheres as the naive counting would 
suggest. We find that the above picture, in which the 
finite frequency vibrational modes are continuously turned 
on as the distortion is increased, unifies the scenario for 
aspherical particles with the one for spheres. 

Since we are ignoring in the analysis below the trivial 
rotations about their symmetry axis, our spheroids actu- 
ally have five rather than six nontrivial degrees of freedom. 
In this analysis, the critical contact number Zc is therefore 
10 rather than 12 [17]. 



Preparation of packings of ellipsoids. In our nu- 
merical study, we simulate ellipsoids that are spheroidal: 
they have two principal axes a and b that are the same 
and a third one c that is different, a = b =/= c. De- 
pending on the ratio of the axes, one distinguishes be- 
tween oblate ellipsoids with aspect ratio e = c/a < 1 
("M&M"s) and prolate ellipsoids with e = c/a > 1 
("cigars"). Particles are of equal size and mass m = I 
and interact with the repulsive Gay-Berne potential [24]: 
V{rij,aij) ^ k/a{{<Tij - rij)/ao)", where = [r^- - r,,] 
is the distance between the centers of the ellipsoids i and 
_j, and aij is the orientation-dependent range parameter. 
This expression is the same as the power law potential 
which has been used in jamming studies of frictionless 
spheres [1]: when Vij = aij the ellipsoids just touch and 
for Vij < aij they repel with a force which is a power of the 
effective overlap aij — r^j; they do not interact otherwise. 
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The range parameter (Jij is defined as 
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Here, u is a unit vector along the principal axis of the 
ellipsoid, li,; = sin 9i cos Lpi x -f- sin 9i sin i^iy -\- cos 9i z, and 
X is the dimensionless parameter, X — (^^ ~ + !)• 

Spatial scales are measured in units of ~ 2a = 2b. For 
anisotropic particles, the parameter k that appears in the 
expression for the potential is known as the well-depth 
anisotropy function and is typically taken to be a function 
of the three directions, k = fc(ui,Uj,fy). In order to 
simplify the expression for the potential, we take k = 1 
(for a one-sided harmonic potential this means that the 
bond stiffness equals 1 for every contact). 

As in previous studies [1], we have used two types of 
interactions, harmonic, with a = 2, and Hertzian, with 
a — 5/2. We focus here on the results for the harmonic 
potential; results for Hertzian forces, are qualitatively the 
same, but will be presented in [25]. The configurations 
were made by initially placing N particles at random in 
a cubic box with periodic boundary conditions. We use 
conjugate-gradient energy minimization to obtain stable 
static configurations, and determine the critical jamming 
density 4>c for each of our initial configurations by in- 
crementally compressing (decompressing) until zero pres- 
sure is reached. We then compress the system to obtain 
zero-temperature compressed configurations. Rattlers are 
removed recursively [25]. We study two system sizes, 
N — 216 and N = 512, and average over about 100 inde- 
pendent initial configurations, at each value of Z and each 
density (/). Here we will show results for N = 216 only. 

From the linearized equations of motion for the coupled 
translational and rotational vibrations of the ellipsoids, we 
construct the dynamical matrix V for each configuration. 
We then diagonalize V, whose eigenvalues are the squared 
frequencies, w^, of the vibrational normal modes. 

Analysis of the spectrum of vibrational- 
rotational modes. — An illustration of how the number 
of rotational, translational and zero-frequency modes per 
particle at jamming varies as a function of coordination 
number Z is shown in Fig. la. There are 3 translational 
modes per particle (red line); in the spherical limit, Z ^ 6, 
there are (10 — 6)/2 ~ 2 zero energy rotational modes 
per particle, because the isostatic number for spheroids 
is 10. As Z increases, the number of non-zero frequency 
rotational modes per particle increases as (.^ — 6)/2 (blue 
line), while the number of zero modes decreases by the 
same amount (gray line). 

In Fig. lb we show the coordination number versus as- 
pect ratio e of spheroids. The black symbols correspond to 
configurations evaluated very close to the jamming thresh- 
old 0c (e) for each value of the aspect ratio, e. The other 



colors correspond to compressions Scj) = (j) ~ 4'c relative 
to the threshold jamming density (pds)- Note that Z de- 
pends both on e and 5(j). The horizontal dashed line at 
(Z — 6) = 4 corresponds to Zc = 10 which is the Maxwell 
criterion for rigidity of spheroids. We have checked in all 
cases that the number of zero-frequency modes per par- 
ticle at threshold is precisely {Z^ — Z)/2 in accord with 
Fig. la; this is shown by the gray crosses. The inset shows 
that for both oblate and prolate spheroids at the thresh- 
old 5Z = {Z -Q) = 6.6(3)|&|°-^°(''), where fc ee e - 1, in 
agreement with results for two-dimensional ellipses [20,26] . 

In Fig. Ic and Id, we show the averaged density of 
states D{Ld) for six typical situations. In Fig. Ic, we show 
D{lo) for spheroids that are close to spheres, Se — —0.04, 
for three different compressions: close to jamming at 
5(f) < 10~^ (black line), at 50 = 10~^ (green line) and 
for relatively large compression, at Scf) = 10~^ (blue line). 
We find that for small 6e, the system behaves nearly as 
if it was made from spheres but with a new "rotational" 
band of excitations. The plateau in the translational band 
of D{uj) still exists with a sharp onset, w*, determined by 
6Z. Our systems are too small to see the elastic plane- 
waves below (jj* . As we will show, uj* scales in the same 
way as the plateau onset for spherical systems [2-5]. The 
rotational band lies below the translational band and ex- 
tends over the range w* < w < < w*, as marked in 
Fig. Ic. As we will quantify, the spectrum is therefore de- 
scribed as having a lower-frequency rotational band sep- 
arated by a gap from w = as well as by a gap from a 
higher-frequency translational band. 

In Fig. Id, we show D{uj) for highly non-spherical parti- 
cles, 6e = —0.33 for the same three values of compression 
as shown in Fig. Ic. For these systems, the gap between 
the two bands has disappeared. Generally, the gap near 
a; = in packings which have large e and Z k, Z^ =10 
opens up with compression above the jamming threshold. 
We illustrate this in the inset of Fig. Id. This is in com- 
plete agreement with the scenario for spheres [1,2]. 

We can determine the nature of the excitations in 
the two bands by analyzing the eigenvectors of the dy- 
namical matrix. We first look at the relative contribu- 
tion of the rotational degrees of freedom to the mode, 
u^{i), where fi = 1,2,3 labels the translations and 

= 4, 5 labels the two Euler coordinates of the orien- 
tation of each particle. In Fig. 2, we plot the rotational 
contribution (ul) = E^i E^=4 "^(O/ E»=i 
(solid line) and the translational contribution (u^) = 

(dashed line) sepa- 
rately. The lower band, existing below lUs, is predomi- 
nantly rotational in nature while the upper band, above 
oj*. is translational. This is most pronounced when e is 
small as shown in Fig. 2a. In the limit as e approaches 0, 
we find that the contribution of (u^) in the upper band 
falls oflt as (3.57(1) • 10"'') w-2 07(i) up to the onset of 
localized modes at high frequencies. The scaling ^ 
is precisely what one expects from perturbation theory if 
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Fig. 1: The average contact number Z and density of states, 
D{lo), of harmonic packings, (a) Illustration of ho-w the num- 
ber of different modes per particle (excluding rattlers) in the 
packings at jamming varies as a function of the average contact 
number Z. Figs, (b-d) show that, as the ellipticity is increased, 
the average contact number increases. For ^ < 9, there are two 
well-defined bands: a rotational band with (Z — 6)/2 modes per 
particle and a translational band with three modes per inter- 
acting particle as is the case for spheres. Upon increasing Z, 
the number of zero modes decreases as zero-modes are con- 
verted into finite-frequency rotational modes. Above Z « 9, 
there is only one band, (b) Z as a, function of the ellipticity 
Se = e — 1 and distance from jamming S(j> for our 216-particle 
packings. The sharp decrease around the spherical case Se — 
is consistent with earlier results [17,21,22] for hard ellipsoids 
and spherically capped rods. The log-log plot of 5Z versus \Se\ 
in the inset shows that the rise of Z at jamming is consistent 
with a 6Z ~ scaling. The crosses in the main plot for 

small values of Se show that twice the measured number of 
zero-frequency eigenmodes per particle plus Z — 6 add up pre- 
cisely to 4, in accord with panel (a), (c) The density of states 
for slightly oblate ellipsoids, Se = —0.04, for our packings close 
to jamming and at two compressions. The existence of two 
bands separated by a gap as well as a gap at zero frequency is 
clearly visible, (d) For larger ellipticities the two bands merge, 
as is illustrated here for Se = —0.33. The inset shows the data 
for Se = 1 where Z ^ Zc = W at jamming (see panel a). In 
accord with this, the gap near zero frequency increases with 
increasing compression (and therefore increasing Z). This is 
consistent with the argument that lj* increases as Z increases 
above Zc- 



the rotational degrees of freedom are weakly coupled to 
the translational ones. Fig. 2b shows the data at a larger 
£ where the bands have just merged. In this case most of 
the modes are of mixed character, as the rotational and 
translational contributions are comparable. 

We can also determine the homogeneity of the modes 
in space by computing the participation ratio P{lo) — 

(Eti E^=i K(')y/NEli EU °f ^^^h mode. 

Here wc will concentrate only on the lower, rotational 
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Fig. 2: (a) Plot of the rotational component {u^) (solid line) 
and the translational component (u^) (dashed line) of the 
eigenmodes for Se — —0.02 as a function of uj. The lower fre- 
quency band is predominantly rotational ((u^) — 1), while the 
upper frequency band is essentially translational {{ut) — 1). 
The red line indicates that at high frequencies rotational con- 
tribution decreases as uj~^. (b) The same as in (a), but with 
Se = —0.17, when the gap between the two bands has just 
closed and most modes have mixed character. 



band, since similar studies of the translational band in 
spherical systems have already been reported [27]. Fig. 3 
shows that at low values of e, the participation ratio is 
small and that for the highest frequencies near lUs the 
modes become highly localized. Due to finite-size ef- 
fects, our present data must be inconclusive; however, we 
surmise that the rotational modes are quasi-localized for 
small 5e in the large system limit. For the two largest val- 
ues of Se shown in this figure, the bands have just merged. 

In Fig. 4a the frequency of the lower edge of the rotation 
band, cu* is plotted vs. \Se\ = |e — 1|. For small |&|, the 
behavior is essentially linear; for large |fe|, when Z at 
jamming approaches 10, the gap closes. 

Fig. 4b shows a;^ and w* as functions of |fc| for har- 
monic configurations of prolate ellipsoids that are close 
to the jamming threshold. We find u* = 1.4(3)|fe|° '^(i) 
and ujs = 3.5(3)|(5e|^-^^^^. The scaling of can be un- 
derstood as the maximum frequency of a libration mode. 
As Fig. 3 shows, this mode is strongly localized, so we 
can obtain the scaling of the maximum frequency by esti- 
mating the torque response for rotating a single ellipsoid, 
keeping the other ones fixed. For a small rotation by an 
angle d9, a contact is compressed or decompressed by an 
amount aQ\6s\d9, where (Tq is the size of the ellipsoids. 
This changes the normal force by an amount kcff(7o\S£\d9, 
with fceff the effective bond strength (for the harmonic 
data fceff = !)• For a slightly oblate or prolate ellipsoid, 
the change in torque is smaller by an amount \Se\, hence of 
order kc{fao\Se\'^d9 as the bond vector and the normal to 
the surface are almost parallel. This implies that the max- 
imum frequency ojs ~ V^eff \5s\. Similar results are found 
in two dimensions [26] . The inset of Fig. 4b shows lUs and 
oj* as functions of |&| for oblate ellipsoids. Values of the 
exponents that we find for oblate ellipsoids are, within the 
error bars, the same as for prolate ones, with prefactors of 
LUs and u!* that are 10% and 15% higher, respectively. 

In Fig. 4c we show uj* as a function of SZ for different 
compressions and ellipticities. u* is still dependent solely 
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Fig. 3: Participation ratio P{uj) of the rotational modes for 
various e as a function of frequency. The eigenmodes at the 
upper edge of the rotational band are seen to be strongly lo- 
calized; throughout the rest of the band -P(w) is quite flat and 
rather small. The data at the largest ellipticities {Se — —0.17, 
and 0.20) correspond to values where the gap between the two 
bands has just closed — the dip in these data is the vestige of 
the merging of the two bands. 



on SZ and therefore the translational band docs not de- 
pend on whether the increase in 6Z occurs due to an in- 
crease in compression or an increase in the aspect ratio of 
the particles. For spheres, the onset of the translational 
band is determined by the excess number of contacts. Our 
results show that for ellipsoids, the same scenario applies, 
irrespective of the origin of the excess contacts. Note that 
the upper and lower limits of the rotational band, cug and 
uj* do not obey this simple behavior but depend differently 
on ellipticity and compression [25]. 

In this study we have dealt exclusively with small sys- 
tems. This raises the question about what finite-size ef- 
fects are important. As mentioned above, one important 
effect is due to plane waves. The lowest-frequency plane 
wave has frequency that is inversely proportional to the 
linear dimension of the system: Lu„iin = CTkmm where 
kmin = 2-K / L with L the linear system size [7]. In the 
systems we have studied, these are predominantly in the 
upper translational band^. In larger systems plane- wave 
elastic modes would begin to populate the gap between the 
rotational and translational band — they would hybridize 
with the modes found in this study, so that localized ro- 
tational modes, for example, would become resonant or 
quasi-localized. 

Conclusions and outlook. — In conclusion, this 
study solves the problem of how the new degrees of free- 
dom associated with non-spherical objects are incorpo- 
rated into the normal-mode spectrum at the jamming 
threshold. Earlier findings [17-20] — that the isostatic 
conjecture breaks down for ellipsoid packings in the regime 
where Z < 10 — suggest that what happens for spheres 
does not immediately apply to more complex shapes. As 
a result, the packing problem of spheres has sometimes 
been viewed as an anomaly [17,23]. If this extended to 
the nature of the jammed state and its dynamic response. 
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^On the basis of Fig. 2 of [7], we estimate that the frequency of 
an elastic mode that occurs at our largest densities A</< is of order 
0.3; since oj* is almost always smaller, in our relatively small systems 
the translational band appears gapped. 
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Fig. 4: Results for the various frequencies defined in the text 
and in Fig. fc. (a) Results for tj^, the lower edge of the rota- 
tional band, as a function of ellipticity at the jamming thresh- 
old. For large ellipticities a;* decreases as Z approaches 10. (b) 
Scaling of ctJ^ and lo* with 5e at the jamming threshold. Data 
is for prolate ellipsoids. lOa , the frequency corresponding to the 
upper edge of the rotational band, which exists for small ellip- 
ticity, scales approximately linearly in Se, in agreement with 
the argument given in the text. Instead lo* scales as j^ej"'^'^-' 
(the red lines are the best fits of the data). The point where 
the two lines cross marks the vanishing of the gap between the 
two bands. Inset: the same data for oblate ellipsoids. Numeri- 
cally, the values are very close to those for prolate ellipsoids at 
the same ellipticity. (c) Scaling of u>* for various compressions 
for oblate and prolate ellipsoids, showing that uj* is determined 
by the contact number only, and that uj* varies linearly in 5Z, 
just as it does for spheres. 



any perturbation from spherical symmetry would quali- 
tatively change the character of marginally-jammed solid 
determined for spheres at Point J. Instead, we find that 
the structure of the normal-mode spectrum remains ro- 
bust. The new modes that are introduced do not affect 
the plateau in the density of states until the spheroid el- 
lipticity becomes large. Moreover, these rotational modes 
appear to be localized so that they should not be efficient 
at transporting heat. The onset of the modes in the trans- 
lational band still depends only on the excess number of 
contacts 6Z = Z — 6 as it does for spheres, irrespective 
of whether the excess contacts result from compression or 
particle asphericity. Thus, the singular jamming transi- 
tion for spheres. Point J, in which the onset of jamming 
coincides with the isostatic point, controls the behavior 
of systems of particles with more complex shapes, just as 
it controls the behavior of sphere packings that are com- 
pressed away from the transition. 

There are two important regimes for ellipsoid packings: 
the first deals with small values of |fe| = \e — 1| where 
the physics is a perturbation around the case of spheres; 
the second deals with large values of |fc|, where Z Zc 
(with Zc = 10 for spheroids and 12 for general triaxial 
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ellipsoids). In the large \Se\ case, the system is well de- 
scribed by a theory in which a plateau in the density of 
states opens up near zero frequency just as it did for the 
case of spheres near Z = 6. An interesting open question 
is how the gap closes when Z ^ 10 for large cUipticitics. 

The insights we obtain from the present study are com- 
plementary to those obtained by including friction in the 
vibrational spectrum of a jammed solid [6]. Here we have 
a situation where the system jams when there are many 
fewer contacts than arc needed according to the Maxwell 
rigidity criterion, while in the case of friction there is al- 
ways an excess of contacts compared to the minimum nec- 
essary for stability [6,14,28-30]. Thus, in the case of fric- 
tion there was never a question of a possible change in the 
underlying picture of jamming threshold. Moreover, at 
the Coulomb threshold to mobilization, the response dic- 
tated by the friction law is inherently discontinuous, which 
makes the properties of packings with friction much more 
sensitive to the preparation history. 

The addition of orientational degrees of freedom does in- 
troduce a new band that is essentially rotational in charac- 
ter while the upper band remains nearly completely trans- 
lational. It is interesting to note that the low-frequency 
rotational modes couple in a simple manner to the higher- 
frequency translations. This is perhaps one of the root 
causes for the very wide spectrum of dielectric response 
that has been observed in glass-forming liquids [31]. The 
rotational modes should appear in the heat capacity. In- 
deed, the boson peak seen ubiquitously in glasses has 
been ascribed [32] to the excess modes associated with 
the plateau of D{ll!) — the fact that the jamming scenario 
is found to be so robust is crucial for its applicability to 
glasses. 
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